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The boundary of a fractionalized topological phase can be gapped by condensing a proper set of 
bosonic quasiparticles. Interestingly, in the presence of a global symmetry, such a boundary can 
have different symmetry transformation properties. Here we present an explicit example of this 
kind, in the double semion state with time reversal symmetry. We find two distinct cases where 
the semionic excitations on the boundary can transform either as time reversal singlets or as time 
reversal (Kramers) doublets, depending on the coherent phase factor of the Bose condensate. The 
existence of these two possibilities are demonstrated using both field theory argument and exactly 
solvable lattice models. Furthermore, we study the domain walls between these two types of gapped 
boundaries and find that the application of time reversal symmetry tunnels a semion between them. 


I. INTRODUCTION 

The interplay of topology and symmetry can lead to 
interesting phenomena in quantum many-body systems. 
In particular, in the presence of global symmetries, one 
topological phases can divide into several different phases 
with the fractional excitations in the system transform¬ 
ing under symmetry in different ways. Much recent ef¬ 
fort has been devoted to the classification of such ‘Sym¬ 
metry Enriched Topological’ (SET) phases by identify¬ 
ing possible ways for the symmetry to act on the frac¬ 
tional excitations 1 11 . One possibility is for the fractional 
excitations to carry fractional quantum numbers of the 
global symmetry. For example, in an electronic system 
composed of charge e electrons, the fractional excitations 
in the v s= 1/3 fractional quantum Hall state can carry 
charge e/3. A systematic counting exists for this class of 
SET phases 3 (although it is not completely clear which 
of these phases exist in pure two dimension and which 
exist as the surface of a three dimensional system). 

Another possibility in SET phases is for the symme¬ 
try to map one type of fractional excitation to another. 
For example, the double semion topological order, which 
exists in for example a double layer fractional quantum 
Hall system with v = ±1/2, is time reversal invariant. 
Time reversal symmetry maps between the semion (with 
topological spin i) and the anti-semion (with topologi¬ 
cal spin —i) while keeping their combination - a bosonic 
quasiparticle - invariant. 

Even though the semion / anti-semion are not individ¬ 
ually invariant under time reversal symmetry, one might 
wonder if we can still make sense of the ‘quantum num¬ 
ber’ or local symmetry transformation on each of them. 
This question becomes more concrete when we consider 
the boundary of the topological state. The boundary 
of the double semion state can be gapped, as shown in 
Fig.l, by condensing the bosonic quasiparticle, which is 
the combination of the semion and the anti-semion. Be¬ 
cause both the semion and the anti-semion obtain non¬ 
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FIG. 1. A double semion state whose boundary is gapped by 
condensing the bosonic quasiparticle. 

trivial phase factors (—1) when braiding around the bo¬ 
son, they become confined in the condensate. Moreover, 
because the bosonic quasiparticle is condensed, i.e. they 
can appear and disappear freely, a semion becomes in¬ 
distinguishable from an anti-semion. Now it becomes 
reasonable to ask about the local time reversal symme¬ 
try transformation of the semion (anti-semion) and there 
could be two options: T 2 = 1 and T 2 = —1. Such a dis¬ 
tinction leads to measurable physical effect in the system. 
Imagine creating a pair of semions in the bulk, separating 
them from each other and bringing them to the gapped 
boundary. When the semion transforms as T 2 = -1, a 
local Kramer degeneracy appears at the location of each 
semion while with T 2 = 1 semion no local degeneracy is 
expected. 

So what does this ‘time reversal quantum number’ 
of the semion imply? Does it label different bulk SET 
phases or does it correspond to different boundary con¬ 
ditions for the same bulk phase? In this paper we show 
that the latter is true. In particular, we use both field the¬ 
ory arguments and exactly solvable models to show that 
there are two types of Bose condensates in the double 
semion state, with different coherent phase factors of the 
condensed boson. In one of them the semion transforms 
as 7~ 2 = 1 while in the other it transforms as T 2 = —1. 

The paper is organized as follows: In section II, 


we present a simple field theory argument for the re¬ 
sult, which is supported by exactly solvable models con¬ 
structed in section III. In section IV, we consider the 
situation where different segments of the boundary of 
the double semion state is gapped in the two different 
ways and ask what happens at the domain wall between 
the segments. We find that the domain wall carries ex¬ 
tra degeneracy protected by time reversal symmetry and 
the symmetry action tunnels a semion between pairs of 
domain walls. In section V, we discuss how this is all re¬ 
lated to the symmetry protected topological phase with 
Z 2 and time reversal symmetry, which becomes the dou¬ 
ble semion state under study by gauging the Z 2 symme¬ 
try. Finally, we conclude in section VI and compare this 
example with similar models studied previously. 


II. FIELD THEORY ANALYSIS 

The double semion state contains an abelian topolog¬ 
ical order with three types of fractional excitations: the 
semion s, the anti-semion s' and their combination - the 
boson b = ss'. The topological spins of the three are 
i, —i and 1 respectively and the mutual statistics be¬ 
tween s and s' is trivial. In field theory language, the 
double semion topological order can be described as a 
U( 1) x U( 1) Chern-Simons theory: 

L = ^-e^aid^al - a\d^al ( 1 ) 

whose edge state can described as 

L e = j^d x (j)id t (j)i - zfcdtfa (2) 

where only the topological term in the Lagrangian is 
shown. 

Time reversal symmetry action on the edge fields 0i 
and <p 2 can be written in two ways: 

7i : 0i -A 02, 02 0i (3) 

Or equivalently, 

72 : 01 02, 02 01 + 7T (4) 

7i and T 2 differ by a gauge transformation 12 

9 • 01 01 + TTj 02 02 (5) 

g acts as g = (—l) Ns where N s is the number of semions 
on the edge and the action is trivial on all local operators 
of the form e z2n( f>i+' l2rn( f >2 w ith integer n and m. 

It appears that a semion, generated by e 1 ^ 1 (or the 
anti-semion generated by e^ 2 ) transforms as T 2 = 1 un¬ 
der 71 and as T 2 = —1 under 72- So which time re¬ 
versal transformation should we use? That depends on 
the boundary condition we choose for the edge theory. 
In particular, the edge state described by Eq.2 can be 


gapped out by adding a Higgs term 

A L = —Acos(20i — 202 + ol) (6) 

with A > 0. When A is large enough, the bosonic quasi¬ 
particle, generated by , is condensed on the edge. 

However, there are two types of condensates which pre¬ 
serve time reversal symmetry, one with a = 0 and the 
other with a = 7 r. 

When a = 0, the term A L has two classical minima: 

0i — 02 = 0, 0i — 02 = 7T (7) 

The two minima are related by g and hence are physically 
identical (no local observable distinguishes them). Each 
minimum is invariant under 7i. Therefore this is the form 
of time reversal transformation that we should consider 
and the semion e 1 ^ 1 (or the anti-semion e 1 ^ 2 ) transforms 
as time reversal singlets T 2 = 1 on the boundary. 

When a = 7 r, the term A L has two classical minima 
as well: 

01-02=7r/2, 01 02 = 7r/2 (8) 

The two minima are related by g again but neither of 
them is invariant under 7i. Instead they are preserved by 
72- Therefore, T 2 is the manifest time reversal operation 
on the edge when boson is condensed with a m 7 r and the 
semion e 1 ^ 1 (or the anti-semion e^ 2 ) transforms as time 
reversal doublets T 2 = — 1 on the boundary. 

Let us look more carefully at how the time reversal 
transformation of the semion / anti-semion depends on 
the coherent phase factor of the Bose condensate. Note 
that in the condensates, the process of creating or an¬ 
nihilating a boson pair is associated with a phase factor 
of 

< >= e ia _ _j_i ^ 

while the process of boson hopping always comes with a 
phase factor of 1 

< e *(0i(®)-02(x)) e -*(0i(x')-^2(a:')) > = 1 (IQ) 

Therefore, the total wave function of the condensate 
reads 

I Ac) = Y,{e ia ) N Y, 1 * 1 . -.*2 n) (11) 

N Xi,...,X2N 

where the inner sum is over all possible position configu¬ 
rations of 2 N bosons and the outer sum is over all integer 
N. Note that there are always an even number of bosons 
because it is a quasiparticle (self) boson and can only be 
created in pairs. 

The relation between the condensed phase and the time 
reversal transformation of the semions can be understood 
as follows: when applying time reversal to a semion, it is 
mapped to an anti-semion and hence a boson is created; 
when applying time reversal again and mapping the anti- 
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semion back to a semion, another boson is created. That 
is, the process of applying T 2 to a semion is accompanied 
by the creation of a boson pair and hence 

T 2 = e ia (12) 

on each semion. 


III. EXACTLY SOLVABLE MODEL 
CONSTRUCTION 

Using the exactly solvable model of double semion in¬ 
troduced by Ref. 13? , we can construct the two types of 
condensates corresponding to a = 0 and a = 7 r and show 
explicitly that the semion excitations transform respec¬ 
tively as T 2 = 1 and T 2 = — 1 under time reversal. 



FIG. 2. Boundary (thick black line) of the double semion 
state (lower half plane) gapped by a = 0 Bose condensate 
(upper half plane) with the semion excitation transforming 
as T 2 — 1 under time reversal. 


A. o = 0 condensate and T 2 = 1 semion 


Consider the double semion model on Honeycomb lat¬ 
tice, with one spin 1/2 living on each link. To simplify 
notation, we write spin operators cr x , a y , a z as A, Y, Z. 
The Hamiltonian in the topological bulk contains vertex 
terms A v and plaquette terms B p 13 . 


H t — — A v + B p 

V p 

where 

a = ru„ z> 


(13) 


(14) 


i G v labels links attached to a vertex v and j G p labels 
links around a plaquette p. rij — 0 if Zj = 1 and nj = 1 
XZj = -l. 

Time reversal T acts as complex conjugation in the 
Z basis. Then the Hamiltonian is time reversal invariant 
and so is the ground state. Note that the form of the pla¬ 
quette term B p used here is slightly different from that 
used in Ref. 13, but they are equivalent when all vertex 
terms A v are satisfied. The form used here is simpler for 
our purpose because it is explicitly time reversal invari¬ 
ant. 

It is easy to check that 1. all Ay s commute, all Ay s 
commute with all B p s and B p s commute with each other 
when the A v constrains are satisfied 2. ( B p ) 2 = 1. 

Now let’s condense the bosonic quasiparticle in the 
double semion state and gap out the boundary. In the 
simplest form, the boson condensation can be achieved 
by enforcing a Zp term on the link variables. That is 


H C0 = -J2 Z k-'E A v ( 15 ) 

k v 

Obviously the second term is redundant. We include it 
here just for comparison with later cases. In the dou¬ 


ble semion state, the Z& term creates / annihilates boson 
pairs and also hops bosons around. Therefore, when the 
Hco term dominates over the Ht term, the boson is con¬ 
densed. Here we consider the exactly solvable situation 
as shown in Fig. 2 where the Hamiltonian in the region 
below the thick black boundary line is Ht and that above 
the boundary line is Hco . All Hamiltonian terms com¬ 
mute with each other. 

From simple counting, we see that the boundary be¬ 
tween the condensate and the topological region is to¬ 
tally gapped. Moreover, this condensate does not break 
time reversal symmetry. A boson-boson pair is created 
by string operator 


Wl =l[z k ( 16 ) 

keL' 

with L' being a string in the dual lattice. W\ obviously 
has eigenvalue 1 everywhere in the condensate. Boson 
hopping is also generated with this term and also comes 
with a phase factor of 1. Therefore, the condensate gen¬ 
erated with Hco is the a = 0 condensate discussed in the 
previous section. 

Now let’s see what happens when we create a semion on 
the boundary. A semion-semion pair is created, as shown 
in Fig. 2 along the red dotted line L. We have chosen a 
particular direction for this string operator. The string 
operator acts as 13 

W s = 

keL Vi v r 

where cq and a r are phase factors in the Z basis at ver¬ 
tices where the string turns left (vi) and right (v r ) re¬ 
spectively. cq = ±1 is always real. a r = (i) Ur acts on 
the leg to the left side of L at this vertex. Obviously in 
the condensate, the string operator costs linear energy 
and the semion is confined. 

Under time reversal symmetry, this string operator 
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changes. In particular, the a r phase factors change to 
ay 1 and the difference is 

wt =i[z k ( 18 ) 

keL' 

where L' is the string in the dual lattice to the left side 
of L, as shown in Eq.2 with the dotted blue line. This 
is nothing but the string operator for creating a pair of 
bosons on the boundary, which has eigenvalue 1. There¬ 
fore, the state with a pair of semions on the boundary is 
invariant under time reversal and each semion transforms 
as a time reversal singlet with T 2 = 1. 


B. a = 7T condensate and T 2 = — 1 semion 



FIG. 3. Boundary (thick black line) of the double semion 
state (lower half plane) gapped by a = n Bose condensate 
(upper half plane) with the semion excitation transforming 
as T 2 = — 1 under time reversal. The Hamiltonian terms in 
the condensate B Pa Z a bB Ph are represented as two red dots 
( B Pa and B Ph ) connected by a grey dashed line via a yellow 
triangle (Z ab ). 

Now let’s define a second type of condensate and see 
how the semion on the boundary can transform as a time 
reversal doublet. The Hamiltonian realizing the second 
type of condensate is 

Hen = - ^ B Pa Z ab B Ph - ^ A v (19) 

<ab> v 

where B Pa and B Pb are neighboring plaquette operators 
sharing a link ab and the sum is over all such pairs. In the 
subspace where all vertex constraints A v are satisfied, it 
is easy to see that all terms in Hen commute with each 
other ( Z ab anti-commutes with B Pa or B Pb ). Moreover 

(■ B Pa Z ab B Ph ) = 1 (20) 

The relative signs in Eq.19 are fixed by the requirement 
that the Hamiltonian is not frustrated. With our choice 


of signs it is not, since 

1 = A v = ( B Pl Zi2B P2 )(B P2 Z23B P3 )(B P3 ZsiB Pl ) (21) 

when both B Pa Z ab B Ph = 1 and A v = 1. 

The B p .ZijB p . terms create boson pairs or hop bosons 
around in this condensate. When this term acts on a pair 
of plaquettes with either 2 or 0 bosons ( BplBp2 = 1), 
there is an extra sign relative to its action on a pair of pla¬ 
quettes with only one boson (. BplBp2 = —1). Counting 
the phase factors at both ends, pair creation / annihila¬ 
tion is accompanied by a phase factor of —1 while boson 
hopping has a phase factor of +1. Note that this is true 
even when we put the small pieces of boson string oper¬ 
ators together and make longer strings. 

(B P1 Z 12 B P2 ) (B P2 Z 2S B P3 )... (B Pm _ i Z(m—l)mB Prri ) 

= B Pl Zi2...Z( rn -i) rn B Prn 

Hence in the ground state we have 

Bp 1 Zi2—Z( m _ 1 ) m B Prn = 1 (23) 

Therefore, the condensate generated with Hen corre¬ 
sponds to the a = 7 r condensate discussed in section II. 
Even though an extra phase factor is present in the wave 
function of this condensate, the condensate is still time 
reversal invariant. 

Now let’s create a semion pair on the boundary along 
the red dotted line in Fig. 3 using the same string opera¬ 
tor W s as in the previous section. Similar to the previous 
case, the string operator W s violates the vertex A v terms 
at its end points, but this has extra consequences in this 
new type of condensate. In particular, when A v is vio¬ 
lated, the plaquette operators B p around this vertex no 
longer commute with each other. To restore exact solv¬ 
ability, certain terms need to be removed, introducing 
local degeneracies into the low energy Hilbert space. As 
we will show below, this local degeneracy indicates the 
presence of a local Kramers pair with T 2 = — 1 under 
time reversal. 





FIG. 4. Semion excitation on the boundary with a = tt con¬ 
densate transforms as a Kramer doublet with T 2 = —1. 

Let us zoom in on the semion as shown in Fig. 4. To 
expose the two fold Kramer degeneracy related to the 
semion, we redefine the Hamiltonian locally as follows: 
imagine breaking the link between plaquette 1 and 5 into 
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two parts and adding a link s sticking into plaquette 5 as 
shown in Fig. 4. The state of the link s is initially set to 
be n s = 0. We can choose the semion string operator as 
going into the condensate and ending on link s. The state 
of the link s is hence flipped to n s = 1. By doing so, we 
have moved the vertex violation to the end of link s which 
does not affect plaquette 1 and 2 but only plaquette 5. 

Due to the existence of link s, B P5 needs to be rede¬ 
fined. According to the string-net rule given in Ref. 13, 
B P5 can be obtained by merging a semion loop into pla¬ 
quette 5. Now with the link s occupied by a semion 
string, direct calculation shows that 

4=*n x » n n njni+i z t (24) 

kep5 j^P5 

where t is the lower half of the link between plaquette 
1 and 5. We can explicitly check that B Ph is Hermitian, 
( B P5 ) 2 = /, it commutes with all other plaquette terms 
already in the Hamiltonian (in the sector where all A v 
constraints are satisfied), but is NOT time reversal in¬ 
variant. Indeed, we find 


Because 

(TQT~ 1 )Q = B Pl Z t B Pl Z t = -1 (28) 

we see that the degenerate states form a local Kramer 
pair under time reversal. Or in other words, on the 
boundary with a = tt condensate, semions transform as 
Kramer doublets under time reversal symmetry. 

IV. DOMAIN WALL BETWEEN TWO TYPES 
OF BOUNDARIES 

Although the two types of condensates, with a = 0 and 
a = 7r, give rise to different boundaries with the double 
semion state, they do not correspond to different phases. 
Indeed, both a = 0 and a = tt condensates are short 
range entangled states with time reversal symmetry. As 
we know that there is no nontrivial time reversal sym¬ 
metry protected topological phase in 2D 14 , both conden¬ 
sates belong to the same phase. Therefore, the interface 
between these two types of condensates can be gapped 
out without breaking time reversal. 


TB^T- 1 = - B P5 (25) 

To indicate the time reversal violation, we denote this 
term with an empty circle in Fig. 4. 

To allow the semion string to end without violating 
time-reversal symmetry, we are therefore forced to re¬ 
move terms containing Bp$ from the Hamiltonian. This 
suggests that four terms need to be removed: B Pi Zi$B P5 , 
i = 1,2, 3,4. In fact, we can recombine the Hamiltonian 
terms and add some terms back (see Fig. 4) 

^ B Pi Zi5Z( i+1 ^B Pi+1 (26) 

i= 1,2,3 

These terms are time reversal symmetric, commute with 
all other terms in the original Hamiltonian and take 
eigenvalue 1 in the ground state. Now counting the num¬ 
ber of terms we find that we are missing one term and 
hence have a local two-fold degeneracy. 

Is this the local Kramer degeneracy we are looking for? 
Obviously this two fold degeneracy corresponds to eigen¬ 
states of B P5 (with eigenvalue ±1) which commutes with 
all other terms in the Hamiltonian but is also indepen¬ 
dent of them. Notice that B P5 anti-commutes with T, 
hence T interchanges the two states with B P5 = =bl. 
Therefore, to determine the T 2 value on this local de¬ 
generacy, we need to find a local operator Q which maps 
between these two states and check its transformation 
under time reversal 2 . That is, Q needs to commute with 
all terms in the Hamiltonian but anti-commute with B P5 . 
One possible choice is 


Q — B Pl Z t 


(27) 



FIG. 5. Domain wall (red dot) between the two types of 
boundary between the double semion state and the a = 0 
and a = tt condensates respectively. 

Now we can ask the question of what happens on the 
domain wall between the two types of boundaries. Imag¬ 
ine a situation as shown in Fig. 5 where a 2D sphere is par¬ 
titioned into three parts, occupied by the double semion 
state, a = 0 condensate and a = tt condensate respec¬ 
tively. The interface between any two parts is gapped 
with time reversal symmetry being preserved and now 
we can investigate the property of the two domain walls 
(red dots in Fig. 5). In this section, we are going to see 
whether there are degeneracies associated with the do¬ 
main walls, and how they transform under time reversal 
symmetry. 


A. Field theory analysis 

From the field theory analysis in section II, we see that 
the boundary with the a = 0 condensate is in state | A) = 
| 0 i — 02 = 0) (or equivalently | A ) = | 0 i — 02 = tt)) 
and the boundary with the a = tt condensate is in state 
| B) = |0i — 02 = 7r/2) (or equivalently \B) = |0i — 02 = 
—7 t/ 2)). The difference between | A) and | A) (or | B) and 
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| B)) is an artificial one as they are related by the gauge 
transformation 


9 • 01 4>1 + 7T, 02 02 

(29) 

In reality, the gauge symmetry is not broken so we need to 
restore the symmetry and write the two boundary states 
as 

\A) + \A) 

(30) 

and 


\B) + \B) 

(31) 


Therefore, the domain walls between the two types of 
boundaries carry a two fold degeneracy and time reversal 
action in this degenerate subspace is equivalent to the 
tunneling of semions between the two domain walls. 

More generally, if the boundary is divided into 2 N al¬ 
ternating segments with 2 N domain walls in between, 
similar analysis shows that there is a 2 2N ~ 1 fold degen¬ 
eracy protected by time reversal symmetry. 

We are going to confirm this conclusion with exactly 
solvable models in the next section. 


B. Interface between the two types of condensates 


When the boundary contains both parts, there are two 
possible configurations, 


|V>i) = | AB) + 1 AB) 

(32) 

|V> 2 ) = \AB) + \AB) 

(33) 


which form a two fold degeneracy on the boundary as 
long as time reversal symmetry is not broken. This can 
be shown as follows. Suppose the two domain walls be¬ 
tween type A and type B boundaries are at L and — L 
respectively, then the local operator near L 

Q 1 = e -i(4> 1 (L+e)-(t>2(L+c)) ^ 4 ) 

(e is small and finite) tunnels a boson across the domain 
wall at L and takes ±1 eigenvalues in the two states. 

OilV’i) = l^i), Oi|V> 2 ) = - 1 ^ 2 ) ( 35 ) 

However, 0\ cannot be added to the Hamiltonian be¬ 
cause it breaks time reversal symmetry, either in the form 
7! or ?2 • 


Tf^OiTi = -Oi, r 2 - 1 0 1 r 2 = 0\ (36) 

On the other hand, the operator 0 2 which tunnels a 
semion from — L to L 


First we need to show that the interface between the 
two types of condensates can be gapped. This is ex¬ 
pected because the two condensates are both time rever¬ 
sal invariant short range entangled states. Because there 
are no nontrivial symmetry protected topological phases 
with time reversal symmetry in two dimensions, the two 
condensates are in the same phase and should be able to 
connect to each other in a gapped way without breaking 
the symmetry 

Consider an interface between the two condensates as 
shown in Fig. 6 



'O. Jt ^ 



FIG. 6. Interface between condensate I and II. 


0 2 = e -iM-L) e iML) ( 37 ) 

maps between |?/q) and ^ 2 ) and anti-commutes with 0\. 
0 2 is a nonlocal operator and cannot be added to the 
Hamiltonian to split the degeneracy. Moreover, as 0\ 
and O 2 generate the full operator algebra of the two di¬ 
mensional space spanned by \^\) and we see that 
|Vq) and 1 ^ 2 ) are necessarily degenerate if time reversal 
symmetry is preserved. 

We can take either form of the time reversal action, 7i 
or 72 , and we find that their action on these two states 
is the same as 0 2 


The upper half plane is in the a = 1 r condensate with 
Hamiltonian terms B Va Z a \ ) B Vh (red dot - yellow square 
- red dot) and the lower half plane is in the a = 0 con¬ 
densate with Hamiltonian terms Z & (yellow square). The 
green squares label the degrees of freedom on the inter¬ 
face and form the interface Hilbert space. 

Let’s be more explicit about this. First, notice that 
the a = 7r condensate can be mapped to the a = 0 
condensate by unitary transformation 


u= n 


(/ + iB p ) 
1 + i 


(39) 


T |Vh) = ^2) 


(38) where the product is over all placates in the condensate. 
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We can see this explicitly from 

(I-\-iB Pa ) (I-\-iBp b ) y (. I—%Bp a ) (I—iB Pb ) 

1-H l+i a6 1-i 1 —i (40) 

= Bp a Z a bB Ph 

Therefore, after applying (7, all links inside the two con¬ 
densates (yellow squares) are in the state | rii = 0) and 
U II |0) is the ground state of these regions. 

Next, we enforce the closed loop constraints 

A v = Y[Zi (41) 

iEv 

along the interface (thick black line) as indicated by the 
green triangles. 

Now we are left with one 2 dimensional degree of free¬ 
dom per plaquette, labelled by k s = 0,1 as shown in 
Fig. 6. \k s = 0) and \k s = 1) are eigenvalue 1 and — 1 
eigenstates of Z s respectively. Each state in the interface 
Hilbert space can be then be written as 

I ^{kg}) = U \ks)\k a+ i)\k a + k 8+1 ) ni°> ( 42 > 

V C 

where the first product is over all downward pointing ver¬ 
tices on the interface and the second product is over all 
the links inside the two condensates (the ones with the 
yellow squares). Note that the links on the interface are 
not all independent due to the A v constraints. There is 
one free degree of freedom per plaquette on the interface. 
As given in Eq.42, each state \^{k s }) is the eigenvalue 
e l7vks eigenstate of iZ s B Ps , where s labels links along the 
interface and p s labels the plaquette on the a = 7 r con¬ 
densate side of link s. Under time reversal, which acts 
as complex conjugation, iZ s B Ps is mapped to —iZ s B Ps . 
Therefore, each \^{k s }) is not time reversal invariant. 

To find a time reversal symmetric interface, we must 
add some terms to the boundary that mix the \^pk s ) 
states, respect T, and gap the interface out. While it 
may be hard to directly find such an operator, we can ap¬ 
ply the unitary transformation U = n P ec„ ApH and 
map the interface Hilbert space to that spanned by 

Km) = II \ k s)\h+i)\k s + k s+ 1) n |0) (43) 

V c 

The transformed interface Hilbert space now takes a sim¬ 
ple direct product form of local degrees of freedom la¬ 
beled by k s and allows simpler analysis of possible Hamil¬ 
tonian terms. Even though the unitary U is not local, it 
preserves the spectrum and hence a gapped edge in this 
basis is also gapped in the original basis. 

There is one complication though: The action of time 
reversal is also transformed under U. Before we can write 
down time reversal invariant Hamiltonians, we need to 
find the correct time reversal transformation in this new 


basis. 



Therefore, the effective time reversal action T on the 
interface Hilbert space is complex conjugation in the 
\ip® k y) basis and 

n*«n(-)"* (i ""' +i) 

Although it looks complicated and non-onsite, we know 
that it should have a short range entangled ground state. 
Indeed we find that under local unitaries 

V = JJ(i)"* (1 “ n * +l) • (46) 

The effective time reversal action is mapped to 

VTV^ = Y[X s K (47) 

which has a gapped symmetric Hamiltonian X s . Note 
that in Fig. 6 X s acts on four links, the two links labeled 
by k s and also the ones labeled by k s -i+k s and fc s +fc s+ i. 
This is illustrated in Fig. 6 with red dotted lines. 

Therefore, the transformed interface Hilbert space 
spanned by | ^ kg y) can be gapped in a time reversal in¬ 
variant way by effective Hamiltonians 

E w ( 48 ) 

In the original interface Hilbert space spanned by |'0{/c s }) 5 
the Hamiltonian reads 

^UV ] X S VU ] (49) 

s 

In appendix A, we explicitly confirm that the term 
UV^X S VW is indeed real. 


C. Tri-junction between a — 0, a — tt condensates 
and double semion 

Now we have found a gapped interface between the 
a = 0 and the a = it condensates, we can study the 
domain wall between the two types of boundaries as the 
tri-junction between the two condensates and the topo¬ 
logical state (double semion). 

As shown in Fig. 7, topological region is on the outside 
and the two condensates are in the middle. Following pre¬ 
vious discussion, Hamiltonian in the topological region is 
given by (red dots and green triangles) 

H t = Y j A v + Y j b p (50) 

V p 
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D z J^A V • B P 

FIG. 7. Tri-junction between a = 0, a = tt condensates and 
the double semion state. 


Hamiltonian in the a = 0 condensate is given by (yellow 
squares) 

h Co = -J2 z ^-J 2 a - ( 51 ) 

k v 

and that in the a = tt condensate is given by (blue dot - 
yellow square - blue dot) 

Hen = - T B Pa Z ab B Pb - A v (52) 

<ab> v 

Moreover, we enforce the closed loop constraints A v along 
the interface between the two condensates (thick black 
line). 

Now we can see what the low energy Hilbert space is 
composed of. Apply again the unitary transformation 
U = ripecq ~ • Note that U commutes with all the 
terms in Ht and the A v constraints on the interface, 
maps the B Pa Z a i ) B Ph terms in Hen to and leaves Hco 
untouched. Therefore, states in the low energy Hilbert 
space can be written as 

n I t+^ p i4fe<>) ( 53 ) 

peCu 

where 

\4>° {ki} ) = I M n \ki)\ki + i)\h + h+ 1) n |0) (54) 

V C 

The first product is over all downward pointing triangles 
on the interface between the two condensates and the 
second product is over all links in the two condensates. 
Note that the two k\ s near the two tri-junctions are the 
same, due to the closed loop constraints over the whole 
state. \(j>Ds) is the wave function in the topological re¬ 
gion (including all the links around red dots as shown 


in Fig. 7). The exact form of | (fros) depends on whether 
ki = 0 or k\ = 1. The two can be mapped into each 
other by running a semion string operator from one tri¬ 
junction to another along, for example, the orange line 
shown in Fig. 7. 

Following a similar line of reasoning as discussed in the 
previous section, we find that the interface can be gapped 
with time reversal invariant terms 

UV ] X S VU ] (55) 

where U is again a product over plaquettes in the a = tt 
condensate U = Y[ peC ^ and V = 

There is one major difference though from the situa¬ 
tion discussed in the previous section. X 1 is now a non¬ 
local operator. Not only does X\ change k\ to 1 — £q, 
it also changes the form of \(/)ds} by running a semion 
string from one tri-junction to another. Therefore, if we 
require time reversal symmetry and locality, we are left 
with a two fold degeneracy with k\ = 0 or 1 respectively. 
Time reversal symmetry maps k\ to 1 — k\ which can 
be equivalently accomplished by X\ in this degenerate 
Hilbert space. Therefore, time reversal symmetry acts 
as semion tunneling in the low energy Hilbert space of 
tri-junctions. 


V. RELATION TO Z 2 x Z^ SYMMETRY 
PROTECTED TOPOLOGICAL PHASE 

The discussions in the previous sections tell us that 
there is one symmetry enriched topological phase for dou¬ 
ble semion topological order with time reversal symme¬ 
try. However, there can be two different gapped bound¬ 
aries with time reversal symmetry. The double semion 
topological order is a twisted Z 2 gauge theory and by un¬ 
gauging the Z 2 symmetry we can get Z 2 x Zj symmetry 
protected topological phases. So how is our conclusion 
about double semion SET consistent with what we know 
about Z 2 x Zj SPT? 

From the classification of SPT we know that there are 
four phases with Z 2 x Zj symmetry 14 . In particular, 
there are two root phases: phase (1,0) which has non¬ 
trivial Z 2 SPT order and trivial action of time reversal 
and phase (0,1) with trivial Z 2 SPT order and projective 
action of time reversal on Z 2 twist defects 15 . Phase (1,1) 
is their combination with both nontrivial Z 2 SPT order 
and projective Z 2 twist defects under time reversal. 

After gauging the Z 2 symmetry, we get SET phases. 
Phase (0,0) and (0,1) lead to the usual non-twisted Z 2 
gauge theory (toric code order) and the resulting states 
differ in the way Z 2 fluxes transform under time reversal 
(as singlet or doublet). Therefore, after gauging, phase 
(0, 0) and (0,1) results in two different SET phases. 

On the other hand, phase (1,0) and (1,1) gauge into 
the twisted Z 2 gauge theory (double semion order). As 
there is only one time reversal SET with double semion 
order, the two SPT phases have to merge into one upon 












gauging. To understand how this happens we notice that 
phase ( 1 , 0 ) and ( 1 , 1 ) can be mapped into each other 
by relabeling the group element. In particular, in the 
symmetry group Z 2 x = {I,9,T,gT}, there are two 
anti-unitary operators T and gT- If we exchange the 
two, the group is still a^x Zj group and phase ( 1 , 0 ) is 
mapped into phase (1,1) and vice verse. Therefore, upon 
gauging the Z 2 symmetry, phase ( 1 , 0 ) and ( 1 , 1 ) become 
the same SET. 

More explicitly, the edge state of the Z 2 x ZJ SPTs 
can be described in general with 

Le = (t>2 (56) 

Different phases correspond to different ways symmetry 
act on 0i and 0 2 . In both phase (1,0) and (1,1), Z 2 acts 
as 

g '• 4> 1 0i + tt, 02 —>■ 02 + 7 T (57) 

Time reversal acts respectively as 

71 : 01 01,02 —^ —02 /ro\ 

72 : 01 01 + 7T, 0 2 ->► -02 

e *Oi +</> 2 )/2 ( or e z( 0 i- 0 2 )/ 2 ^ crea ^ es a twist defect on 
the edge and transforms as Tj 2 = 1 and 7^ 2 = — 1 . 

Combining g with 7i we get 

gTl : 01 -► 01 + 7 T, 02 -02 - 7 T (59) 

Redefine 0 2 = 0 2 + 7t/ 2, we find 

^7l : 01 01 + 7T, 02 -02 (60) 

which is exactly the same as the action of 7^ • Hence after 
gauging 7i and 72 are equivalent, and both correspond 
to the same SET phase. However a memory of these 
two distinct SPT phases is retained at the boundary of 
our SET, where depending on boundary conditions the 
semion can transform with either T 2 = 1 or — 1. 


VI. CONCLUSION AND DISCUSSION 

To summarize, in this paper we have learned the fol¬ 
lowing things about the double semion topological order: 

1. There are two different ways to condense boson 
and gap out the boundary while preserving time 
reversal symmetry in the double semion state. One 
corresponds to a Bose condensate with a coher¬ 
ent phase factor 1 and the semion excitations on 
the boundary transform as T 2 = 1. The other 
corresponds to a Bose condensate with a coherent 
phase factor of —1 and the semion excitations on 
the boundary transform as T 2 = — 1. 

2. A pair of domain walls between the two types of 


boundaries carry a two fold degeneracy. Time re¬ 
versal symmetry acts in this two dimensional space 
by tunneling a semion from one domain wall to an¬ 
other. 

3. There is only one SET phase with time rever¬ 
sal symmetry and double semion topological or¬ 
der. The different transformation properties of the 
semions under time reversal symmetry is a pure 
boundary effect. In the bulk, semions are mapped 
to anti-semions and it is not meaningful to talk 
about the time reversal representation carried by 
semion itself. 

4. Different SPT phases can become the same SET 
phase after partly gauging the unitary symmetry 
of the system. Examples of this kind have been 
pointed out in Ref. 16. 

These results can be generalized to other SET or SPT 
phases. In appendix B, we present another example of 
SPTs coalescing upon gauging with unitary symmetry. 
We leave the study of more general cases to the future. 

We want to comment briefly on the relation between 
this double semion example and some previous studies of 
gapped boundaries of topological states. 

A simple yet very interesting case was the boundary 
state of toric code topological order. Even in the ab¬ 
sence of symmetry, there are two types of boundaries 
corresponding to the two types of (self) bosons in the 
toric code 17 . Each domain wall between the two types 
of boundaries carry a Majorana mode, giving rise to a 
2 n fold degeneracy for N + 1 pairs of domain walls (with 
fixed fermion parity). 

Similarly, degeneracies arise with domain walls in our 
double semion example. With N pairs of domain walls, 
there is a 2 2N ~ 1 fold degeneracy. Of course, this degener¬ 
acy requires the protection of time reversal symmetry and 
can be completely removed by adding time reversal sym¬ 
metry breaking local terms. If time reversal symmetry 
is preserved, the degenerate states can only be mapped 
to each other through non-local operators which tunnel 
semionic excitations from one domain wall to another. 

The existence of such domain wall degeneracy has 
also been noticed in (fractional) topological insula¬ 
tors where domain walls between ferromagnetic gapped 
edges and superconducting gapped edges carry Majorana 
(parafermion) zero modes 18 22 . Such models are differ¬ 
ent from the double semion example studied here in that 
symmetry is broken in order to gap out the edge. 

Finally, Wang and Levin studied different ways to gap 
out the edge of a ‘strong pairing insulator’ 23 , either with 
an interface with a topological insulator or with an in¬ 
terface with a trivial insulator. Symmetry is preserved 
in these two kinds of edges, but the interface between 
the two condensates - the topological insulator and the 
trivial insulator - is always gapless when symmetry is 
preserved. 
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FIG. 8 . Each X s term flips four links a, /3, 7 and S. 


where 


S a = (A2) 


and 


V* 


•na-i+l-n^i r\n a -in a + (l-np)(l-np +1 ) 


(A3) 


Complex conjugating UV^X S VW then gives 


(UV^XsVW)* 

_ 1 — iBpi 1 — iB P2 1 +iB P2 

- i +r~ 

= UB Pl B P 2 a*<T*a%a* s V;B Pl B P2 W 


(A4) 


where p\ and P 2 are the plaquettes above link a — 1, a 
and /?,/?-hi respectively, as shown in Fig. 8. 

Therefore, UV^X S VU^ is real if 


O TD _CC X X X T/* D TD — X_.X_X_.XT T 

J^Pi G OL a /3 a 8 V(f) £>Pi Up2 — °cx a f3 a 5 

This can be explicitly checked as 

b p , b P2 n v x = n^(-r“-+^ +i 
(-) n “^ +l V ; = v, p 

and B P1 B P2 commutes with V^. Therefore, 


(A5) 


(A6) 


(uv^x s vu ] y = uv ] x s vu ] (A7) 

Appendix A: Explicit calculation of T-invariance of 
the gapping terms between the two condensates 


In this section, we are going to explicitly verify that the 
UV^X S VU^ term given in section IV B is indeed real. 

First, let’s write the X s term in full as cr^cr^cr^cr|, as 
shown in Fig. 8 Then 

VA.V 

= St t _ 1 Slo%oZ<7%<7tS a - 1 S f) (Al) 

= ( 7 * 0^0 


Appendix B: Example of SPTs coalescing upon 
gauging with unitary symmetry 

Consider SPTs with x (unitary) group in 2d. 
The classification is Z\. The root phases have the edge 
structure, 

L = ^-d x ^d T e (Bl) 


10 








and transformation properties: 


Phase (1,0,0) 


.v(i) . 

.z 2 . 

v( 2 ) . 

z 2 . 


cj) — y cf> 7T, 6 — y 6 -\- 7r, 

trivial 


Phase (0,1, 0) :Z^ : 

v( 2 ) . 


trivial 


6^6 + : 


Phase (0, 0,1) 


.z 2 . 

v( 2 ) . 

z 2 . 


cj) — y (/> 7T, $ —^ $ 

0 —>> 0 , + 7r 


Combining phases (0, 0,1) and (0,1,0) we obtain 
Phase (0,1,1) : Z^ : 0-^0 + 7 r, 6^6 

zf> ■. <t> -s- 4> + 7T, e ->■ e + tt(B2) 


Indeed, denote the edge modes of (0, 0,1) as 0i, #i and 
edge-modes of (0,1, 0) as 02 , # 2 - Then the combined edge 
is described by, 

£ = (B3) 


actions of #2 and #3 interchanged. Now imagine gauging 
the Z^ group in the two phases ( 0 , 0 , 1 ) and ( 0 , 1 , 1 ). 
The action of gi in both cases is identical and the twist 
defect is given by e l9 = e l9 ! 2 . In terms of 6 the edge 
theory is, 


L = ^d x cbd T e (B 6 ) 

Z7T 

i.e. after gauging we get a toric code topological order. 

( 2 ) 

Now, Z\ ; remains a global symmetry in the resulting 
SET, 

( 0 , 0 , 1 ) with Z^ gauged; Z 2 2 ^ : 0 —>► 0 , ^ —>► (9 -|- 7 t /2 
(0,1,1) with Z^ gauged; Z 2 2 ^ : 0 —)> 0 + 7 r, + 7 t /2 

(B7) 

Clearly, as SETs the two phases are the same since all 
local degrees of freedom transform in the same way (the 
transformation properties differ by a pure “gauge” trans¬ 
formation 0 0 + 7 r). 

We can also think about gauging the remaining Z 2 ; 
symmetry. In the case of (0, 0,1) the twist defect of g<± is 
e uf) _ e ^/2 Thus, we get an overall Z 4 topological order, 


Let’s define 0' x = 0i, 6 [ = 61 — 62 , 0 2 = 0 1 + 02 , 6' 2 = 62 . 
The action has the same form in the primed variables as 
in the unprimed. The transformation properties of the 
primed variables are, 

: 0i 0i + tt, 6[ —y 6[ 

02 02 + n j 62 —y 62 

: 0 i 0 i j #1 “^ 

02 “^ 02 + ^ 5 #2 “^ $2 + ^ 

(B4) 

Adding a term —A cos( 6 ^) we gap out the 0 ' l5 ^ modes. 
The transformation properties of 02 , #2 are then exactly 
the same as in Eq. (B2).] 

Let us denote the generator of Z 2 1 ^ as c/i, the generator 
of Z^ as g 2 and g% = # 1^2 • Then, in phase (0,1,1) under 
#3 we have 

gs : 0 —)> 0, 6 ^ 6 + it (B5) 

So the phase ( 0 ,1,1) is like the phase ( 0 , 0 , 1 ) but with the 


L = (B 8 ) 

Z7T 

where e 1 ^ is the twist defect of Z 2 2 ^ and e l9 is the twist 
defect of Z^\ 

Now, in the case of (0,1,1) the twist defect of $2 is 
e 1 ^/ 2 e l9 / 2 = e^e 10 . If we are interested in the “overall” 
topological order, we can still use the e 1 ^, e l ° basis, ob¬ 
taining a Z 4 topological order. However, the way that 
the twist defect of $2 is “embedded” within this topo¬ 
logical order is different. In (0,0,1) it is just (which 
is a boson), whereas in ( 0 , 1 , 1 ) it is e^e 19 (which is a 
semion). If we are treating the two phases as SET phases 
(with Z 2 1 ^ “fully gauged”), then the difference of an ex¬ 
tra factor of e l9 is irrelevant, since the twist defect of the 
global symmetry Z\ ; can always trap an extra anyon 
(vison) of the SET, e l °. Thus, we cannot distinguish the 
two SET phases. However, if we are working with SPT 
phases and putting the twist defects in “by hand”, then 
we know whether the twist defect of gi is present or not, 
so can distinguish the two SPT phases. 
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